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Abstract 

The cross-section ratio = a(e + e~ — ► bb)/cr(e + e~ — > qq) and the bottom and charm forward- 
backward asymmetries Ap B and Ap B are measured using event samples collected by the OPAL 
detector at centre-of-mass energies between 130 and 189 GeV. Events with bottom quark pro- 
duction are selected with a secondary vertex tag, and a hemisphere charge algorithm is used 
to extract Ap B . In addition, the bottom and charm asymmetries are measured using leptons 
from semileptonic decays of heavy hadrons and pions from D* + — > D°7r + decays. The results 
are in agreement with the Standard Model predictions. 
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1 Introduction 



Many measurements of heavy quark production have been performed in e + e~ collisions on the Z 
resonance |1[ . Among the measured parameters are the production cross-sections of bottom and 
charm quark pairs relative to the hadronic cross-section, and R c , and the forward-backward 
asymmetries A^ B and Ap B . In this paper, measurements of R^, Ap B , and Ap B at energies above 
the Z° resonance are presented for e + e~ — > Z°/ / y* — > qq events, where the effective centre-of- 
mass energy \fs ! after initial-state radiation is required to satisfy a/ s' / s > 0.85. Similar 
measurements have been performed previously by other collaborations |2|,[3|]. 

The data collected with the OPAL detector at LEP at centre-of-mass energies between 
130 GeV and 189 GeV are analysed. The basic techniques are similar to those adopted in 
previous OPAL measurements The R^ measurement is based on the selection of a sam- 

ple enriched in bb events obtained with a secondary vertex tagging technique. Both rip and 
rz information^ from the silicon microvertex detector are used in the selection of bb events, 
improving on |4], [|] where only rip information was used. This provides higher efficiency at 
comparable purity, thus enhancing the statistical precision of the measurement. 

Two measurements of forward-backward asymmetries are performed. In the events selected 
by secondary vertex tagging, the asymmetry Ap B is measured using a hemisphere charge tech- 
nique to identify the direction of emission of the primary quark. An independent measurement 
of both the bottom and charm forward-backward asymmetries is performed using leptons from 
semileptonic decays of heavy hadrons and pions from D* + — > D°7r + decays. The asymme- 
tries obtained with the hemisphere charge method and the lepton and slow pion methods are 
combined. 

In Section ^, a brief description of the OPAL detector and the event selection is given. 
The i?b measurement is discussed in Section ^ followed by a description of the asymmetry 
measurements in Section |J Systematic errors on both the i?b and asymmetry measurements 
are given in Section ||, and in Section || all results are summarised. 

2 The OPAL Detector and Event Selection 

A detailed description of the OPAL detector can be found elsewhere JFj . For this analysis, the 
most relevant parts of the detector are the silicon microvertex detector, the tracking chambers, 
the electromagnetic and hadronic calorimeters, and the muon chambers. The microvertex 
detector is essential for the reconstruction of secondary vertices. The central detector provides 
precise measurements of the momenta of charged particles by the curvature of their trajectories 
in a magnetic field of 0.435 T. In addition, it allows an identification of charged particles through 
a combination of the measurement of the specific energy loss dE/dx and the momentum. The 
electromagnetic calorimeter consists of approximately 12000 lead glass blocks, which completely 
cover the azimuthal range up to polar angles of | cos#| < 0.98. Nearly the entire detector is 
surrounded with four layers of muon chambers, after approximately one metre of iron from the 
magnet return yoke, which is instrumented as a hadron calorimeter. 

1 The OPAL coordinate system is defined as a right-handed Cartesian coordinate system, with the x axis 
pointing in the plane of the LEP collider towards the centre of the ring, the z axis in the direction of the 
outgoing electrons, and 8 and <p defined as the usual spherical polar coordinates. 
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Starting in 1995, the LEP experiments have collected data at increasing energies well above 
the Z° peak. In this paper, the energy points are classified in five different sets, at centre-of- 
mass energies which will be generically called a/s=133, 161, 172, 183, and 189 GeV. Table |1] 
shows the luminosity-weighted mean centre-of-mass energies at which data were taken, and 
the corresponding integrated luminosities. Additionally, calibration data taken at the Z° peak 
during 1996, 1997, and 1998 are used to cross-check the analyses. 



Dataset 


133 GeV 


161 GeV 


172 GeV 


183 GeV 


189 GeV 


(Vi> (GeV) 

Integrated luminosity (pb 1 ) 


133.3 
9.9 


161.3 
9.6 


172.1 
9.9 


182.7 
55.8 


188.6 
177 



Table 1: The effective centre-of-mass energies and integrated luminosities of data collected 
above the Z° peak. The uncertainties on these quantities have a negligible effect on the analyses 
presented in this paper. 



Hadronic events, e + e~ — > qq, are selected based on the number of reconstructed charged 
tracks and the energy deposited in the calorimeters. The selection of the subsample of non- 
radiative hadronic events, defined by the requirement \J s' / s > 0.85, and the identification and 
rejection of W + W~ background are described in detail in ||. The remaining contamination 
from radiative hadronic events with true effective centre-of-mass energy below 0. 85y/s is 5 — 10%, 
depending on the centre-of-mass energy. The residual contamination from four-fermion events 
(mainly W and Z pairs) is largest at yfs = 189 GeV, where it is about 8%. These backgrounds 
to non-radiative hadronic qq events are accounted for in the measurements. 

Jets are reconstructed using the JADE algorithm || with the E0 recombination scheme 0, 
keeping the invariant mass cut-off x m i n = 49 GeV 2 fixed at all centre-of-mass energies. 

Hadronic events are simulated using the PYTHIA Monte Carlo generator [ 10| . Heavy quark 



fragmentation is modelled according to the scheme by Peterson et al. |IT[ with fragmentation 
parameters tuned according to the results in |12[ . Four-fermion background events are simulated 
with the grc4f generator ||13|| . The events are passed through a detailed simulation of the OPAL 



detector 14] before being analysed using the same procedure as for the data. 



3 Measurement of 

The tagging of bb events is based on the long lifetime (~1.5 ps) and hard fragmentation 
of b-flavoured hadrons, which give rise to secondary vertices significantly displaced from the 
primary vertex. The secondary vertex tag described in Section [JTT] allows a clean and efficient 
reconstruction of bb events. In Section the measurement of with vertex tagged events 
is described. 



3.1 Secondary Vertex Tag 

The algorithm used for secondary vertex reconstruction is described in [pjj. For the results 
presented here, a three-dimensional vertex tagging algorithm is used, which takes advantage of 
the precise z information provided by the OPAL microvertex detector. 
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The primary vertex in each event is reconstructed as described in ||16|| , incorporating the 
average beam spot position determined from the measured tracks and the LEP beam-orbit 
measurements as a constraint. Although the beam spot is less precisely determined at energies 
above the Z° resonance than at the Z° peak, the resulting error on the primary vertex position 
is still small compared to the error on the reconstructed secondary vertex position. The angular 
acceptance is restricted to | cos6^| < 0.9, where 6t denotes the polar angle of the thrust axis 
of the event. Charged tracks used to reconstruct secondary vertices are selected as described 
in [|j] and least three of them are required to form a vertex. 

Each hadronic event is divided into two hemispheres by the plane perpendicular to the 
thrust axis and containing the nominal interaction point. For each reconstructed secondary 
vertex the signed decay length L is defined as the distance between the secondary and the 
primary vertex. L is taken to be positive if the secondary vertex is in the hemisphere pointed 
at by the momentum vector of the jet which contains the vertex, and negative otherwise. 
The decay length significance L/ai is defined as the ratio of the decay length and its error. 
Secondary vertices with L/a L > 8 are used to tag bb events. This cut represents the best 
compromise between tagging efficiency and purity. 



3.2 Determination of from Vert ex- Tagged Events 

For the measurement of R^, the number of events tagged by a secondary vertex is determined 
and corrected for tagging efficiency and background. Due to the limited statistics compared 
with the data collected at the Z° peak, a double tag technique as e.g. in cannot be ap- 
plied. In order to reduce the sensitivity of the analysis to the detector resolution, a folded tag 



technique |16| is used: a hemisphere is assigned a tag if it contains a secondary vertex with a 
decay length significance L/<jl > 8, or an anti-tagQ if it contains a vertex with a decay length 
significance L/<jl < —8. The number of anti-tagged hemispheres is then subtracted from the 
number of tagged hemispheres. After subtraction of the four-fermion background, the differ- 
ence between the number of tagged and anti-tagged hemispheres, N — N, in a sample of iVhad 
hadronic events can be expressed as 

N-N = 27V had [(e b - e h )R h + (e c - e c )R c + (e uds - e uds )(l - R h - R c )) , (1) 

where (eb — Cb), (e c — e c ), and (e u d s — e u d s ) are the differences between the tagging and anti- 
tagging efficiencies for a given quark flavour. The tagging efficiencies for uu, dd, and ss events 



are averaged (uds), since they are very similar. For i? c , the prediction of ZFITTER |18[ is used. 
The four-fermion background is determined from Standard Model production cross-sections 
and from selection efficiencies estimated from Monte Carlo simulation. The b-purity for the 
sample with the folded tag is defined as the fraction of bb event hemispheres contained in the 
sample N — N, and is about 75%. The distribution of the decay length significance L/aL is 
shown in Figure [l] for events at 189 GeV centre-of-mass energy, together with the expectation 
from the Monte Carlo simulation. The relative difference between data and Monte Carlo is 
below 1% for the folded-tag rate and (8 ± 4)% for the anti-tag rate, which is more sensitive 
to modelling of the detector resolution. For Lj<j L < —10, the number of events predicted by 
Monte Carlo differs by two standard deviations from the number of events in the data. This 



2 In the literature, tagged and anti-tagged events are sometimes referred to as "forward" and "backward" 
tagged events. This convention is not used here to avoid confusion with the distinction between forward and 
backward event hemispheres in the asymmetry analyses. 
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indicates an incomplete simulation of the detector resolution, which is taken into account in 



the systematic errors as described in Section |5.1.6| below. The agreement between data and 
simulation is similar at the other centre-of-mass energies. 

The selected event sample contains a 5 — 10% contamination of radiative hadronic events, 
with one or more energetic photons emitted in the initial state. In these events the effective 
centre-of-mass energy is reduced to values below a/ s' / s = 0.85 where the predicted value of 
i?b is U P to 30% larger. In addition, the selection of non-radiative events is around 3% less 
efficient for bb final states than for other flavours, because of a larger missing energy due to 
neutrinos in semileptonic b hadron decays. The results are corrected for these effects, which 
are estimated from Monte Carlo simulation. The measured values of are also corrected for 
interference between initial and final-state radiation as described in 0. This correction comes 
from the fact that Monte Carlo, which is used to model the data, does not contain interference 
between initial and final-state radiation. 

The numbers of selected events and of tagged and anti-tagged hemispheres are given in 
Table [2|. The differences in hemisphere tagging efficiency, also listed in Table [|, have been 
determined from Monte Carlo simulation. Their errors include all the systematic uncertainties 



that will be described in Section |5.1| . No systematic uncertainties other than those due to Monte 
Carlo statistics and detector resolution are assigned to the efficiencies in uu, dd, and ss events, 
as they represent a small fraction of the tagged sample. The systematic error is dominated by 
the uncertainties from the event selection, the modelling of b and c fragmentation and decay, 
and from the simulation of the detector resolution. 



Energy 


Events 


N 


N 


Tagging efficiency 
differences 




pSM 
n b 


133 GeV 


745 


153 


11 


e' h = 0.347 ±0.026 
< = 0.045 ±0.009 
e' uds = 0.0035 ±0.0014 


0.190±0.023±0.007 


0.184 


161 GeV 


347 


58 


4 


e' h = 0.331 ±0.013 
< = 0.045 ±0.002 
e' uds = 0.0090 ±0.0010 


0.195±0.035±0.007 


0.171 


172 GeV 


228 


23 


3 


e' h = 0.338 ±0.012 
e' c = 0.041 ±0.002 
e ' uds = 0.0066 ±0.0010 


0.091±0.034±0.005 


0.168 


183 GeV 


1186 


232 


25 


4 = 0.351 ±0.009 
e' c = 0.048 ±0.002 
e ; ds = 0.0084 ±0.0008 


0.213±0.020±0.009 


0.165 


189 GeV 


3209 


551 


87 


e' h = 0.363 ±0.008 
e' c = 0.053 ±0.002 
e' uds = 0.0099 ±0.0008 


0.158±0.012±0.007 


0.164 



Table 2: The numbers of selected non-radiative events, tagged (N) and anti-tagged (N) 
hemispheres, and tagging efficiency differences e' q = (e q — e q ) in the i?b analysis. The errors 
on the tagging efficiencies include systematic uncertainties. The i?b results with the statistical 
and systematic errors and their Standard Model expectations are given in columns 6 and 7. 
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Figure 1: The decay length significance distribution for selected events at 189 GeV centre-of- 
mass energy. Only the most significant value is used for each event. The points with error 
bars represent the data. The dark-shaded histogram represents the expected contribution from 
four-fermion background, and the light-shaded area corresponds to the background from non- 
bb two-fermion events. The unshaded histogram indicates the bb-content in the sample. The 
arrows show the position of the cut. 
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The dependence of the result on the assumed value of R c can be parametrised as 

AR h = b (R c - i? c SM ) . (2) 

The parameter b has been determined separately for each centre-of-mass energy. Its values are 
given in Table |3|. 



Dataset 


133 GeV 


161 GeV 


172 GeV 


183 GeV 


189 GeV 


b = (AR h ) / (R c - R^) 


-0.12 


-0.11 


-0.11 


-0.12 


-0.13 


pSM 


0.223 


0.244 


0.249 


0.253 


0.255 



Table 3: The dependence of the measured R^ value on the assumed value of R c for each 
centre-of-mass energy. The parameter b gives the change in R^ for an assumed deviation of R c 
from its Standard Model prediction. 

As a cross-check, the analysis is repeated on calibration data collected at the Z° peak. A 
value of Rb(\/s = m z o) = 0.221 ± 0.002(stat.) ± 0.010(syst.) is obtained. The systematic error 
has been determined as for the high-energy data samples. This result agrees within the errors 
with the LEP1 combined value of R® = 0.21664 ±0.00076 |L| and can be regarded as a check of 
the evaluation of systematic errors for the measurements at energies above the Z° peak. Note 
that the 4% systematic error which is assigned to the i?b measurement at a/s = 189 GeV is 
larger than the 2% difference with respect to the LEP1 average which is observed at = m z o. 



4 Measurement of Forward- Backward Asymmetries 

For the measurement of heavy quark forward-backward asymmetries, it is necessary to dis- 
tinguish the event hemispheres of the primary quark and antiquark in addition to the quark 
flavour tagging. Two complementary techniques are used. The first analysis provides a mea- 
surement of Ap B for the events that have been tagged by the presence of a secondary vertex. A 
hemisphere charge method is used to distinguish between quark and anti-quark hemispheres. 

The second technique is used for a simultaneous measurement of Ap B and A^ B . It is based 
on the identification of leptons from semileptonic decays of heavy hadrons ( "prompt leptons" ) 
and pions from D* + — > D°7r + decays ("slow pions"). The charge of these particles provides a 
clean distinction between the primary quark and anti-quark hemispheres. 

In Section O , the measurement of A^ B with the hemisphere charge technique is described. 



After a discussion of the lepton and slow pion identification in Section [4.2| and the flavour 
separation of the tagged samples (Section |4.3|) , the asymmetry measurement based on prompt 
leptons and slow pions is presented in Section [4.4| . Finally, the combination of the two mea- 
surements is treated in Section Ol 



4.1 Measurement of Ap B with a Hemisphere Charge Technique 

To obtain a sample enriched in bb events suitable for the Ap B measurement, the secondary 
vertex algorithm used for the i?b analysis is employed. The analysis is limited to the range 
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I cos6*t| < 0.9. Only tagged events are used, since Monte Carlo studies show that with the cur- 
rent statistics, the folded tag would result in a larger error on the forward-backward asymmetry. 
Thus, an event is considered if it contains a secondary vertex with a decay length significance 
L/cjl > 8 and no other secondary vertex with L/<jl < —8. This cut is chosen to minimize the 
overall error. 

Each non-radiative hadronic event is divided into two hemispheres by the plane perpendic- 
ular to the thrust axis that contains the nominal interaction point. Monte Carlo simulation 
shows that the direction of the thrust axis is a good approximation of the direction of emission 
of the initial qq pair. For each hemisphere, the hemisphere charge Qhem is computed as 

f\ _ ^2i=l \Pi\ ' Qi /o\ 
Vhem - j^n ' ^ > 

where the sum runs over all n tracks in the hemisphere, pi is the momentum component of track 
i along the thrust axis, Qi denotes its charge, and k = 0.4 is a parameter tuned on simulated 
events for an optimal charge identification of a primary b or b. To ensure a good hemisphere 
charge reconstruction, only events with more than three charged tracks per hemisphere are 
used. Events are classified according to the sign of the difference Qf~Qb between the forward 
(Qf) and backward (Qb) hemisphere charges, where the forward hemisphere is defined as the 
one that contains the momentum vector of the incoming electron. 

In order to ensure that each event is used at most once in this analysis and the measurement 
based on prompt leptons and slow pions described below, every event is assigned a figure of 
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merit V^ g defined as 



P^' = F h ■ (2P b - 1) , (4) 



where Fb denotes the estimated bb purity as a function of decay length significance and Pb 
stands for the estimated probability of correct charge identification^ as a function of \Qf~ Qb\, 
both determined from the simulation. Events are rejected from the vertex-tagged sample if they 
are also tagged by the presence of a prompt lepton or slow pion with a corresponding figure of 



merit V^g > Vs^, where Vgfg is determined according to Equation [L^ (see Section [O 



Slg Slg 5 " Slg 

below). At y/s = 189 GeV, the events that contain both a tagged secondary vertex and a 
prompt lepton or slow pion correspond to 56% of the vertex tagged and 28% of the lepton 
or slow pion tagged samples, respectively. Of these common events, 58% are assigned to the 
secondary vertex tagged event sample by the procedure described above. It has been checked 
that systematic cross-dependences between the two asymmetry measurements, which may in 
principle be introduced by this method, are negligibly small. Note that the quantity vf&^ is 
not used in the fit which determines Ap B , but only to define the selected sample of events. 

For the final vertex-tagged event sample, the tagging efficiencies for each flavour are deter- 
mined from Monte Carlo and are shown in Table |j. The uncertainties include the systematic 
errors, which will be discussed in detail in Section [5.1| . They are dominated by uncertainties 
in the bottom and charm physics modelling and detector resolution, as in the analysis. 
The b purity is about 60%, with a fraction of four-fermion background up to 5%. The lower 
b purity with respect to the folded tag purity does not limit the precision of the asymmetry 
measurement with the present statistics. Because of the rejection of some of the events that 
are also lepton or slow pion tagged, the final bb efficiencies and purities are lower than those 
which could otherwise be obtained for a vertex-tagged event sample. 

3 The charge identification probability Pb does not depend significantly on the decay length significance of 
the tagged vertex. 
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Flavour q 


Event selection 


Charge identification 


efficiency e q 


probability P q 


d 


0.032 ±0.005 


0.673 ±0.021 


u 


0.036 ±0.005 


0.743 ±0.015 


s 


0.035 ±0.005 


0.683 ±0.020 


c 


0.097 ±0.005 


0.645 ±0.011 


b 


0.383 ±0.025 


0.683 ±0.012 



Table 4: Event tagging efficiencies e q and charge identification probabilities P q at y/s = 
189 GeV, as determined from Monte Carlo simulation. The efficiencies and charge identification 
probabilities are computed for the final selected samples after rejection of events with V^ 3 ^ > 

Vgte . Similar values have been obtained for the other centre-of-mass energies. 



The quantities P q are the probabilities for the hemisphere charge method to correctly iden- 
tify the event hemisphere into which the primary quark q was emitted. They are determined 
from Monte Carlo simulation and are given in Table f| for the different quark flavours at 
y/s = 189 GeV. Similar values are obtained at different centre-of-mass energies. Their er- 
rors include systematic uncertainties which will be discussed in detail in Section |54]. The 
largest contribution to the systematic error of the b and c quark charge identification probabil- 
ities arises from the modelling of heavy flavour fragmentation and decay. For light flavours the 
uncertainties on fragmentation are expected to have a small effect on the total systematic error, 
and only Monte Carlo statistics and detector resolution are considered. Possible detector biases 
in the charge identification probability for positive and negative quarks have been investigated 
using Z° calibration data, and have been found to be negligible. 

The quantity 

x = -sign(Q F - Q B ) ■ | cos6> T | (5) 

is computed for each event, where 9t denotes the polar angle of the thrust axis. Its observed 
distribution at 189 GeV centre-of-mass energy is compared with the Monte Carlo prediction 
in Figure |2|. When only vector and axial vector couplings of the quarks to a gauge boson 
exchanged in the s-channel are allowed, the observed angular distribution of the primary quark 



can be expressed as |i9fl 



— = Ce(x)(l + x 2 + U^x), (6) 

where the quark masses have been neglected. The constant C is for normalization, and e(x) is 
the tagging efficiency as a function of cos 9t for an event. It is assumed that the efficiencies 
are symmetric functions of x, and it has been checked in the simulation that their dependence 
on x is the same for all primary flavours. For other event types (e.g. four-fermion events), the 
predicted differential cross-section is not a second-order polynomial, but the resulting effects 
are negligible within the precision of the measurements presented here. 

Using Equation || the observed asymmetry Ap£ is obtained by maximising the log likelihood 

N N r 8 i 

ln£ = J2 1» {Cefa)} + J2 ^ 1 + + o^fe , (7) 
j=i i=i ^ J 



11 



X 90 

d 

g 80 

d 

70 



t Data 
I I 2-fermions 
J 4-fermions 



OPAL 



60 



50 



40 



30 



20 



10 



-1 -0.7S -0.5 -0.25 0.25 0.5 0.75 1 



Figure 2: The distribution of the observed x = — sign(<5_F — Qb) ' \ cos6 T \ at 189 GeV centre- 
of-mass energy. The points with error bars are the data, and the histogram is the Monte 
Carlo expectation. The dark shaded histogram represents the expected contribution from four- 
fermion background. The histogram showing the Monte Carlo expectation has been scaled to 
the same number of entries as found in the data. 

where the sum is over all N selected events, and A^ B is the only free parameter in the fit. The 
first term is a constant for a given set of events. 

After four-fermion background subtraction, Ap B is determined from the relation 

A&= £ Sq -F q .(2P q -l)-AV (8) 

q=u,d,s,c,b 

Here, s q is —1 (+1) for up-type (down-type) quarks, and A^ B is the forward-backward asymme- 
try for flavour q. The asymmetries for non-b events are fixed to their Standard Model values as 
calculated by ZFITTER. The fractions F q of events of flavour q in the sample are determined 

as 

^=v%r> ( 9 ) 

where R q is the ratio of the cross-section of quark type q to the total hadronic cross-section, 
determined from ZFITTER, and e q is the tagging efficiency determined from Monte Carlo. 
The factor (2P q — 1) is to account for charge misassignment. The small contamination from 
four-fermion events is evaluated as for the i?b measurement and subtracted from the sample. 
Its observed asymmetry is found to be consistent with zero within the available Monte Carlo 
statistics. 

The numbers of tagged events, the observed asymmetries A^ B , and the corrected asymme- 
tries Ap B with their statistical and systematic errors are given in Table |5], together with the 
Standard Model expectations at the different centre-of-mass energies. Most of the systematic 
errors are in common with the i?b analysis, as discussed below in Section |5.1|. The largest 



systematic errors arise from uncertainties in the detector resolution and the event selection 
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procedure. Uncertainties in the fragmentation of light quarks are assumed to be negligible and 
are not considered. For all centre-of-mass energies, the statistical error is dominant. 

While the observed asymmetries are well within the physical range of —0.75 < A^ B < 
0.75, values of Ap B outside this range are possible because of the corrections due to sample 
composition and charge identification probability. All corrected A^ B values are compatible with 
a value inside the physical range. No constraint is applied to force the corrected Ap B values to 
lie within this range, in order to facilitate the combination with the values determined in the 
lepton and slow pion analysis and with measurements by other experiments. 



Energy 


Events 


A obs 
^FB 


^FB 


4 b,SM 
^FB 


133 GeV 


75 


0.12 ±0.12 


0.67 ±0.50 ±0.06 


0.48 


161 GeV 


35 


-0 17 +0 - 18 

U- 1 ' -0.16 


-0.59 toll ±0-04 


0.55 


172 GeV 


14 


0-32 ±gS 


1.9 ±1.3 ±0.1 


0.56 


183 GeV 


157 


0.182 +g;g« 


1.15 toir ±0.08 


0.57 


189 GeV 


372 


0.194 ±g;gsp 


1.22 toil ±0.08 


0.58 



Table 5: The numbers of tagged events with the vertex tag analysis, the observed asymmetries, 
the resulting Ap B values, and their Standard Model predictions. The first error on Ap B is 
statistical, the second systematic. 



4.2 Identification of Leptons and Slow Pions 

Prompt leptons from semileptonic decays of heavy hadrons provide a means of tagging both bb 
and cc events that is largely independent of the secondary vertex tag. In addition, slow pions 
from D* + — > D°7r + decays are used for tagging heavy flavour events. Both prompt leptons and 
slow pions allow a clean identification of the event hemisphere that contains the primary quark. 



4.2.1 Electron Identification 

Electron candidates are required to have a momentum of at least 2 GeV. Tracks with less than 
20 dE/dx samplings in the tracking chamber are rejected to ensure a good measurement of the 
specific energy loss. The difference between the measured energy loss and that expected for an 
electron, divided by the measurement error, is required to be between —2 and 4. 

In this sample of preselected tracks, electrons are identified with the help of an artificial 



neural network, which is described in detail in [I7J]. In addition to the electron preselection, a 



network output A/" c i > 0.9 is required. At 189 GeV, this selection has an efficiency for prompt 
electrons of approximately 25%, defined with respect to the total number of prompt electrons 
that are reconstructed as tracks in the detector. The resulting sample is 75% pure in electrons. 
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After this selection, electrons from photon conversions are an important background in the 
sample. A separate artificial neural network is used to identify pairs of conversion electrons @, 
O ]. The contribution from photon conversions is reduced by requiring a network output of 
A/" cv < 0.4. In the Monte Carlo simulation at 189 GeV centre-of-mass energy, 89% of the 
electrons from photon conversions are rejected by this cut, while 90% of the prompt electrons 
are kept. In Figure [| the A4i and A/" cv output distributions are shown for tracks at 189 GeV 
centre-of-mass energy. 

4.2.2 Muon Identification 

The muon selection proceeds in two steps. First, muon track segments are formed from the 
hits in the muon chambers. Tracks from the central tracking chambers with a momentum 
greater than 2 GeV are extrapolated to the muon chambers. For each track segment in the 
muon chambers, only the "best matching track" is considered for use in the asymmetry fit. It 
is defined as the extrapolated track that has the smallest angular separation a to the muon 
track segment in question. 

In a second step, an artificial neural network trained for muon identification is used to 
enhance the purity of the muon sample. The network uses the following eleven inputs: 

• Information from the matching: 

— The square root of the x 2 for the position match in 9 and <p between the extrapolated 
track and the associated muon track segment in the muon chambers, as described 
in 0; 

— the ratio of distances R mis = a^/at 2 ' of the best and second best matching track 
to the muon segment; this is a measure of how ambiguous the choice of the best 
matching track was in the preselection; 

— the x 2 probability for the matching computed using both position and direction 
information for the track in the central detector and the associated muon track 
segment. 

• Information from the hadron calorimeter: 

— The number of calorimeter layers in the cluster associated with the central track; 

— the number of the outermost such layer; 

— the x 2 probability for the match in 9 and (f) between the track (extrapolated to the 
hadron calorimeter) and the associated cluster. 

• Specific energy loss: 

— The muon dE/dx weight for the track, which is a measure of the probability that 
the track is compatible with a muon hypothesis; 

— (TdE/dx, the error on the dE/dx measurement; 

— the momentum of the track. 

• Geometrical information: 

— The position in | cos 9\ and where the extrapolated track enters the muon chambers. 

The distribution of the neural network output A/^ is shown in Figure |] for "best matching" 
tracks according to the definition above. Muon candidates are retained if A/" M is larger than 
0.65. In Monte Carlo simulated events at 189 GeV centre-of-mass energy, the muon selection 
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Figure 3: The M c \ distribution is shown in figure (a) for tracks at 189 GeV centre-of-mass 
energy that pass the electron momentum and dE/dx cuts. In (b), the J\f cv distribution is given 
for those tracks that also pass the requirement of jV e i > 0.9. In each case, the points with error 
bars represent the data distribution, and the histograms the Monte Carlo simulation, scaled 
to the same number of entries. The dark and light grey areas are the expected contributions 
from prompt electrons and conversion electrons, respectively, while the open area corresponds 
to hadrons. The arrows indicate the accepted regions. 
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results in an efficiency of 43% for prompt muons, defined with respect to all prompt muons 
that are reconstructed as tracks, and a muon purity of 73%. 




0.2 0.4 0.6 0.8 




Figure 4: The A/" M distribution for best matching tracks at 189 GeV centre-of-mass energy that 
pass the muon momentum cut. The points with error bars represent the data distribution, and 
the histogram the Monte Carlo simulation, scaled to the same number of entries. The shaded 
area is the expected contribution from true muons. The accepted region is indicated by the 
arrow. 



4.2.3 Preselection of Slow Pion Candidates 



Pions from D* + — > D°7r + decays, denoted 7r s in the following, are selected based on the kinematic 



properties of this decay. Due to the low momentum, p* = 39 MeV [21], of the decay products 



in the D* + rest frame, pions from this decay have momenta smaller than 

p™* = (E* + p*) = 0.0458 sfe (10) 

2m D *+ 



in the laboratory frame, where E* = (p*) 2 + m?. + , and m^,*+ and m n + denote the D* + and 
7r + masses, respectively. In addition, slow pions have a transverse momentum with respect to 
the D* + flight direction of at most p*, and are thus dominantly found in the core of the jet 
containing the D* + meson. 

Slow pion candidates are required to have a momentum between 1.0 GeV and p™ ax and, 
if at least 20 dE/dx samplings are available, a specific energy loss dE/dx whose probability 
compatibility for the pion hypothesis exceeds 2%. Tracks that form single charged particle jets 
are rejected. The D* + flight direction is estimated by the jet direction, which is recalculated in 
an iterative procedure similar to the one described in [^2|, based on the rapidities of the tracks 



and clusters in the jet containing the slow pion candidate. If the jet mass exceeds 2.3 GeV, the 
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track or calorimeter cluster in the jet with the smallest rapidity with respect to the jet axis is 
removed from the calculation, and the direction is recomputed. The transverse momentum p t 
of the slow pion candidate is calculated with respect to this jet direction. The p\ distribution 
is shown in Figure [5] for tracks at 189 GeV centre-of-mass energy. Slow pion candidates are 
accepted if p\ < 0.02 GeV 2 . 

In Monte Carlo simulated events at 189 GeV, this preselection is 56% efficient and yields 
a sample that contains 5.7% of slow pions. This sample is further enriched with a cut that is 
described in Section 0~3| below. 




0.02 0.04 0.06 0.08 0.1 

p 2 t in GeV 2 

Figure 5: The p 2 distribution for tracks at 189 GeV centre-of-mass energy that pass the slow 
pion momentum and dE/dx cuts. Tracks that form single charged particle jets have been 
excluded. The points with error bars represent the data distribution, and the histogram the 
Monte Carlo simulation, scaled to the same number of entries. The dark grey area is the 
expected contribution from true slow pions from cascade b — > [c,c] — > tc s decays, and the light 
grey area represents c — > tt s decays. The arrow indicates the accepted region. 



4.3 Flavour Separation of the Lepton and Slow Pion Samples 

Three different sources of prompt leptons are considered: b — > £, meaning leptons from semilep- 
tonic decays of b-flavoured hadrons; cascade bottom decays, which include the contributions 
from b —y c — ► £ and b — ► c — > £ processes; and c — > £, leptons from semileptonic decays 
of charm hadrons. The background can be classified as "non-prompt" leptons, i.e. all other 
leptons that are not produced in the decay of b- or c-flavoured hadrons, and particles that are 
mis-identified as electrons or muons. 

For electrons and muons, separate artificial neural networks have been constructed with the 
aim of separating prompt b — > £ decays, cascade b — > [c, c] — > £ decays, prompt c — > £ decays, 
and all other contributions. The technique used is similar to the one described in ||, but the 
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inputs and training have been re-optimized for centre-of-mass energies above the Z° resonance; 
also an artificial neural network for the identification of cascade decays is included. Candidates 
that form single charged particle jets are rejected since they are expected to be dominantly 
produced in leptonic W decays. 

The first network, denoted At>, is designed to separate b — > I decays from all other contri- 
butions. Two networks, Abe and Ac, have been trained on a Monte Carlo simulation that does 
not contain b — > i decays in order to classify lepton candidates that are background to the At> 
net. All three networks use the following input variables, where jet variables are defined with 
the lepton candidate included in the jet: 

• p, the lepton track momentum; 

• pt, the transverse momentum of the lepton track calculated relative to the jet which 
contains the track; 

• L/ul-i the decay length significances of the secondary vertices (if existing) in the jet 
containing the lepton and the most energetic jet in the hemisphere not containing the 
lepton, where secondary vertices are reconstructed with the same algorithm as described 
in Section [TT| 

• the jet charges of the jets containing the lepton and the most energetic jet in the hemi- 
sphere not containing the lepton, each multiplied by the lepton charge, where the jet 
charge is defined as in Equation || with k = 0.4, but using only tracks associated to the 
jet; 

• the forward multiplicity in the lepton jet, defined as the number of tracks with an impact 
parameter significance with respect to the primary vertex larger than 2. For each track, 
the impact parameter is defined as the distance between the primary vertex and the track 
at its the point of closest approach; the impact parameter significance is defined as this 
distance divided by its error; 

• the | cos^l of the jet momentum vector, where 9 is the jet polar angle; 

• the outputs A"ei and M cw of the electron identification network and the conversion finder 
network, respectively, in the case of electrons; 

• the output Aju of the muon identification network, in the case of muons. 

From the Ab, Abe, and Af c network outputs, the following quantities are computed which 
are related to the probabilities of a lepton candidate to come from one of the three sources: 

( b-*[c,c]-^) = n - AH [A~ bc (l -A~ c )] 

1 h) [A bc (l-A'e)] + [A'e(l-A'be)] + [(l-A;)(l-Ab c )] 1 > 



[Ae(l-Ab c )] 



Slg ' ' [Abe (1 " K)] + [K (1 - Abe)] + [(1 " Ac) (1 - Abe)] ' 

The difference in the treatment of the Ab output from Abe and A" c is due to the fact that b — ► t 
decays were omitted in the training of the latter two networks. Only candidates that satisfy 
the condition 

*g = \/(^!r") 2+ ( p r M " , ) 2+ (^") 2 > 01 (12) 

are used in the subsequent analysis. Candidates with lower values of vjfg are expected to be 
dominantly background and to have a negligible contribution to the overall result. Note that 
the quantity "P sig , although used to define the selected sample of events, is not itself used in 
the fit that determines the bottom and charm asymmetries. 
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Similarly to the lepton case, slow pion candidates are classified as cascade b — > [c, c] — > tt s 
decays, c — > ir s , and background. Another two artificial neural networks, JVb c and J\f c , have 
been trained for the separation of the three components in the preselected slow pion sample. 
The following inputs are used by both networks: 

• p, the slow pion track momentum; 

• pf, the transverse momentum squared of the slow pion track calculated relative to the jet 
direction obtained with the same rapidity based algorithm that is used in Section [4.2.3 

• S^-jet, the total energy of the jet containing the slow pion 



E-Ks-sub-jeti the energy of the sub-jet p3j containing the slow pion: Each jet containing a 



slow pion is split into two sub-jets, where the slow pion sub-jet is seeded by the slow pion 
track. In an iterative procedure, any particle that forms a smaller opening angle with 
the slow pion sub-jet than with the remainder of the jet is then assigned to the slow pion 
sub-jet; 

• L/uli the decay length significance of the vertex in the jet containing the slow pion, if a 
vertex is found; 

• the jet charge of the jet containing the slow pion, calculated with k = 0.4, multiplied by 
the slow pion charge; and 

• the | cos#| of the jet momentum vector. 

In addition, the network for identification of b — > [c, c] — > n s decays uses: 

• L, the decay length of the vertex (if existing) in the jet containing the slow pion; 

• the jet charge of the jet containing the slow pion, calculated with a different parameter 
k = 2.0, multiplied by the slow pion charge; 

• the forward multiplicity in the jet containing the slow pion, defined as above; 

Pt)]et-, the scalar sum of the transverse momenta relative to the jet axis of all tracks 
in the jet; and 

• the maximum longitudinal momentum component of any track in the jet containing the 
slow pion, measured relative to the jet direction. 

The quantities V^T^'^ 7 '^ and V^ 77 ^ are computed in analogy to Equation [TT] for prompt 

lepton candidates, while V^^ ns ^ is set to zero. The same requirement on (computed 
according to Equation [12]) is imposed as for prompt leptons. 



If more than one lepton or slow pion candidate per event passes the selection, the one with 
the highest V^ n ^ value is taken. Furthermore, this value is required to be larger than 



the value of v!?^ as determined for the hemisphere charge measurement (see Section [4.1|) if 



the event is also tagged by the presence of a secondary vertex. 

A breakdown of the composition of the samples of electron, muon, and slow pion tagged 
events at 189 GeV centre-of-mass energy together with the efficiencies is given in Table []. In 
Figure D, the output distributions of the flavour separation networks are shown for the events 
that pass all cuts. 



4.4 Measurement of Ap B and Ap B with Leptons and Slow Pions 

The forward-backward asymmetries for bottom and charm, Ap B and A^ B , are extracted from 
the data using an unbinned maximum likelihood fit to the charge signed polar angle distribution 
of the thrust axis T in lepton and slow pion tagged events. It is assumed that T is the axis along 
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Figure 6: The A/b, A/" c , and Abe distributions for lepton and slow pion tagged events at 189 GeV 
centre-of-mass energy. The points with error bars correspond to the data. The open histogram 
represents the Monte Carlo expectation, where the contributions from signal events are shown 
hatched as explained in the figure. The purity in the slow pion selection is limited, which is 
reflected in the fact that the range of M c and Abe values for slow pions does not extend up to 1. 
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Source 


Sam 
Electrons 


Die comp< 
Muons 


)sition 
Slow pions 


Electrons 


3fncienci( 
Muons 


3S 

Slow pions 


b -> [c,c] -> £/7T s 
C -> £/7T s 


19% 
9% 
19% 


19% 
11% 
24% 


3% 
13% 


26% 
16% 
18% 


50% 
31% 
38% 


6% 
13% 


Non-prompt leptons 
Other background 


15% 
37% 


20% 
25% 


85% 





Table 6: The estimated composition of the lepton and slow pion tagged event samples at 
189 GeV centre-of-mass energy after the selection. The efficiencies of the prompt lepton and 
slow pion selection are given for the final event sample with respect to all prompt leptons and 
slow pions that are reconstructed as tracks in the detector. 



which the primary quark-antiquark pair is emitted. The quantity y = — q cos 9t is computed 
event by event, where q is the charge of the lepton or slow pion, and the thrust direction T is 
defined such that T-p > 0, with p being the momentum of the jet containing the lepton or slow 
pion. The inclusive y distributions for electron, muon, and slow pion tagged events at 189 GeV 
are presented in Figure |7|. 

In the fit, both the bottom and charm asymmetries are to be determined. Therefore, the 
fit uses for each event its probabilities to be a correctly tagged bb or cc event or background, 
as determined from the simulation: The events are divided into several subsamples according 
to their At,, M>o N c , and \y\ values. Three bins are used for each quantity, making a total of 
81 subsamples separately for both electron and muon tagged events and 27 subsamples for the 
slow pion tagged events. These subsamples have different bottom and charm purities and are 
fitted simultaneously. 

The cross-section for producing a qq pair is assumed to depend on y according to 

^- 1+ y 2 + l A ^y- (is) 

The likelihood £ su b for one subsample is given by 



E ^U + y 2 + l^y\ , (14) 



hl£sub — 

candidates 

where the sum is taken over all candidates in the subsample. The total likelihood is then given 
by the product of the likelihoods of all subsamples. The expected observed asymmetry Ap B in 
each subsample is computed as 

5 

A^(JV b , jV bc ,jV c , \y\) ^MMhJ^M, \v\) ^fb • (15) 

i=i 

In this equation /, denotes the predicted fraction of leptons or slow pions from source i, and 
A l FB is the corresponding asymmetry: 
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Figure 7: The y distributions of (a) electron, (b) muon, and (c) slow pion tagged events at 
a centre-of-mass energy of 189 GeV are given by the points with error bars. Also plotted is 
the expectation from the simulation together with the contributions from background (open 
histogram) and b — > £ (cross-hatch), b — > [c, c] — > £/ir s (wide diagonal hatch), and c — > 
£/ti s decays (narrow diagonal hatch). Note that in the fit, the tagged events are divided into 
subsamples of different bottom and charm purities to allow an extraction of both the bottom 
and charm asymmetries. 



where x° is the effective B — B mixing parameter. The fractions /, have been calculated from 
the Monte Carlo simulation and depend on the mis-identification probability in each bin of Ab, 
Abe, A/" c , and \y\, on the production rates of bottom and charm quarks, on the semileptonic 
branching ratios of heavy hadrons, and on the hadronisation fractions f(b — > [c,c] — > D* 1 * 1 ) 
and f(c — > D* + ). Variations in the sample composition with |cos#t| are taken into account 
since the fi are binned in \y\. As described in Section |4.1| , this likelihood fit has the advantage 
that the Icos^rl dependence of the efficiency for identifying leptons and slow pions is not 
needed explicitly. The effective mixing parameters x^ ff , x^ ff , and x^ ff are determined from the 
simulation for each selected subsample of events. The mixing parameter \ = 0.118 ± 0.006 
used in the simulation for inclusive bb events is taken as an external input. The small non-zero 
contributions to the observed asymmetry from prompt leptons and slow pions from radiative 
and four-fermion events are accounted for, but left out of the above list for simplicity. The 
assumption that the backgrounds from mis-identified leptons and slow pions do not contribute 
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to the observed asymmetry has been checked and will be discussed in Section 



5.2.2 



The fit is done separately for the data taken at 133 GeV, 161 GeV, 172 GeV, 183 GeV, 
and 189 GeV. In Table [?|, the numbers of lepton and slow pion tagged events are given, and 
the results of the fit for the bottom and charm asymmetries are summarised for each energy 
point together with the errors and correlations. The systematic errors have been evaluated as 



described in Section 5.2 



Energy 


Electrons 


Muons 


Slow 
pions 


Measured 
asymmetries 


Correlation 


Predicted 
asymmetries 


133 GeV 


48 


69 


82 


^B = "0.01 t° H ±0.18 
A C FB = 0.50 t°oil ±0.13 


±14% 


A b F >£ M = 0.48 
A C F - B SM = 0.69 


161 GeV 


16 


37 


42 


4ra= °- 18 -oil ±0.15 
A FB = 0.83 t° H ±0.12 


±14% 


^ B SM = 0.55 
A C FB SM = 0.69 


172 GeV 


14 


29 


20 


A FB = 0.55 ±°:p ±0.18 
A C FB = 0.67 t° f 4 ±0.12 


±23% 


^ B SM = 0.56 
^• B SM = 0.67 


183 GeV 


77 


104 


126 


A FB = 0.59 toll ±0.12 
A FB = 0.56 t° H ±0.11 


±23% 


^ B SM = 0.57 
,4 F < B SM = 0.67 


189 GeV 


190 


269 


359 


A FB = 0.28 ±0.21 ±0.12 
^fb = 0.52 t° ™ ±0.11 


±19% 


<' B SM = 0.58 
A F - B SM = 0.66 



Table 7: For each centre-of-mass energy, the numbers of lepton and slow pion tagged events, 
the results for the bottom and charm asymmetries as measured with the lepton and slow pion 
tag, and their correlation are listed. The first error is statistical, and the second systematic. In 
the last column, the Standard Model predictions for the asymmetries are given. 



As a cross-check, the fit is also performed on the calibration data taken at the Z° peak in 
the years 1996, 1997, and 1998. The results of this fit are consistent with the average of LEP1 
and SLD measurements given in [I]. 

4.5 Combination of the Ap B Measurements 

By construction, the samples used in the hemisphere charge analysis do not have any events 
in common with the lepton or slow pion tagged samples. The bottom asymmetry as measured 
with the vertex tag depends on the assumed value of the charm asymmetry, while the fit to the 
lepton and slow pion tagged events yields results for the bottom and charm asymmetries with 
a non-zero correlation. This dependence on the charm asymmetry has to be taken into account 
when combining the two bottom asymmetry measurements. The vertex tag measurement is 
applied as a constraint in the fit to the lepton and slow pion tagged events that is described in 
Section |4.4j . Since the bottom and charm asymmetries are correlated in the fit to the lepton 
and slow pion tagged events, both the fitted bottom and charm asymmetries are expected to 
change in the constrained fit; but the error on the charm asymmetry is essentially unchanged. 
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The results for A FB and A FB are listed in Table § and are also shown in Figure |8| together with 
the Standard Model expectations as a function of the centre-of-mass energy. Good agreement 
is observed between the measurements and the predictions from ZFITTER. 



Energy 


Measured 
asymmetries 


Correlation 


Predicted 
asymmetries 


133 GeV 


A FB = 0.19 ±0.30 ±0.12 
A FB = 0.50 t°oil ±0.13 


±20% 


4< B SM = .48 
A£ B SM = 0.69 


161 GeV 


A FB = -0.03 t° il ±0.11 
A FB = 0.87 ±{j;*j ±0.12 


±23% 


^ B SM = 0.55 
yl^ B SM = 0.69 


172 GeV 


A FB = 0.82 t° f 2 ±0.14 
A C FB = 0.69 t°oil ±0.12 


±27% 


^ B SM = 0.56 
^| M = 0.67 


183 GeV 


A FB = 0.77 ±g;|t ±0.10 
A C FB = 0.55 ±°;|| ±0.11 


±34% 


^ B SM = 0.57 

y^ M =0.67 


189 GeV 


4hj= 0.63 ig;^ ±0.10 
A FB = 0.50 ±0.11 


±29% 


^ B SM = 0.58 
^ B SM = 0.66 



Table 8: The results of the combined A FB and A FB measurements together with the Standard 
Model predictions from ZFITTER. The first error is statistical, and the second systematic. 



In all asymmetry measurements, R c and Rb are fixed to the Standard Model values as 
predicted by ZFITTER. The dependence of the measured values on the assumed values of i?b 
and R c is parametrised as 



Ai? h . _b/ D x AR r 



AA FB = a\R h )^ + a h (R 



-rtv, it. 



3M" 



A^ B = a c ( J R b )^±a c ( J R c )f^ 



(17) 



where Ai? q = R q — i?q M for q=b, c. The Standard Model values for R h and _R C are given in 
Table |5] together with the values of the coefficients a q (R q >). 

Any systematic error from a common source (see Section |5] for the description of systematic 
uncertainties) is treated as fully correlated between the two measurements. However, there 
are large contributions to the systematic error that affect only one of the analyses, such that 
systematic errors are clearly no limitation to the combination of the two results. 



5 Systematic Errors 

In this section, the evaluation of the systematic errors for the analyses presented here is dis- 
cussed. Both the secondary vertex and the lepton and slow pion tags depend crucially on the 
knowledge of the detector resolution for reconstructing charged particle tracks. Most of the 
other main systematic errors are independent for the two measurements. 
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Figure 8: The R^, Ap B , and Ap B measurements are compared to the Standard Model predic- 
tions. The statistical and the total errors of the measurements are indicated. The measurements 
at y/s ~ m^o have been taken from |), 17, 24, 25]; here, the measurement errors are smaller than 
the points. The behaviour of the curve at y/s ~ 110 GeV is due to the cut y/s'/s > 0.85. 
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energy 




dSM 
ll c 


a b (R h ) 


a h (R c ) 


a c (R h ) 


a c (R c ) 


133 GeV 


0.184 


0.223 


-0.24 


+0.12 


+0.06 


-0.28 


161 GeV 


0.171 


0.244 


-0.07 


-0.06 


+0.14 


-0.86 


172 GeV 


0.168 


0.249 


-0.64 


+2.48 


+0.07 


-0.11 


183 GeV 


0.165 


0.253 


-0.40 


+0.21 


+0.07 


-0.36 


189 GeV 


0.164 


0.255 


-0.38 


+0.23 


+0.10 


-0.30 



Table 9: The Standard Model predictions for the partial hadronic decay widths R^ and R c and 
the dependence of the results on their values are given. The coefficients apply to the combined 
measurements. The symbols are defined in the text. 



5.1 



Systematic Errors on i?b and Ap B with the Secondary Vertex 
Tag 



In Table [K| a breakdown of the systematic error is given for both the R\, and Ap B measurements 
with the vertex tag at 189 GeV centre-of-mass energy. The systematic errors considered for 
these measurements are described in the following paragraphs. 



Error source 


AR h 




Event selection 


0.0038 


0.036 


Final state QCD corrections 




0.011 


b fragmentation 


0.0003 


0.004 


b decay multiplicity 


0.0027 


0.018 


b hadron composition 


0.0017 


0.011 


b lifetime 


<0.0001 


<0.001 


c fragmentation 


<0.0001 


0.001 


c decay multiplicity 


<0.0001 


0.003 


c hadron composition 


0.0011 


0.011 


c lifetime 


<0.0001 


<0.001 


Four-fermion background 


0.0009 


<0.001 


Monte Carlo statistics 


0.0015 


0.031 


Track reconstruction 


0.0051 


0.057 


Total systematic error 


0.0074 


0.078 



Table 10: Systematic error breakdown for the i?b and the Ap B measurements with the vertex 
tag at y/s = 189 GeV. Similar uncertainties have been determined for the other centre-of-mass 
energies. 
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5.1.1 Event Selection 



The bias on the measurement of from the selection of non-radiative events has already 



been discussed in Section ^2. A correction has been applied, and an uncertainty of 100% is 
assigned to the correction. The uncertainty on the corrections for the interference between 
initial and final state radiation results in a systematic error of 0.5% on R\ } . The contamination 
from radiative events has the effect of decreasing the measured asymmetry by 3%. The value is 
corrected accordingly, and a 3% systematic error is assigned. The likelihood fit is based on the 
assumption that the shape of the tagging efficiency as a function of | cos#p| is the same for all 
flavours, which might not be true at the edges of the acceptance. The effect has been estimated 
by dividing the sample in bins of | cos#p|, over which the above assumption is valid, determing 
the asymmetry in each bin independently and comparing their average with the reference result. 
It has been found to affect the asymmetry by less than 1%, and is thus neglected. 

5.1.2 Final State QCD Corrections 

Final state QCD corrections are included in the calculation of the predictions with which the 
asymmetry measurements are compared. However, the experimental event selection is less 
efficient for events with hard gluon radiation due to the cuts on the decay length significance 
(for the vertex tag) and the lepton or slow pion momentum (for the lepton and slow pion 
analysis). The fraction of the QCD correction that has to be applied to the measurements has 
been determined previously at LEP1 and is typically between 0.3 and 0.6 |26|]. The overall 



QCD corrections at a/s = 189 GeV have been determined from ZFITTER to be 0.015 for Ap B 
and 0.022 for Ap B , respectively. Half of this correction is assigned as a systematic error. 

5.1.3 Bottom and Charm Physics Modelling 

Uncertainties in bottom and charm fragmentation and decay properties have been treated as 
follows: 

• b fragmentation: Although the mean scaled energy 2(E^) / ' ^Js of weakly decaying b 
hadrons is expected to change from LEP1 to LEP2 energies, the Peterson fragmentation 
parameter ep, which describes one step in the fragmentation process, is assumed not to 
vary with energy. Simulated bb events are reweighted within the range of 0.0030 < ep < 
0.0048 JT2J, which corresponds to a variation of the mean scaled energy 2(E\ > )/y/s of 
weakly decaying b hadrons in Z° decays in the range of 2{E^) / ^ = 0.702 ± 0.008 fl2| . 

• b decay multiplicity: The charged decay multiplicity of hadrons containing a b quark 
is varied in the Monte Carlo simulation according to |L2] . 

• b hadron composition: The tagging efficiency differs for the various b hadron species. 



The fractions of b hadrons and their errors have been taken from pi |. The fractions 
/(B° + B + ) and /(B°) are varied independently within their errors, and their variation is 
compensated by the b baryon fraction. 

b hadron lifetimes: The lifetimes of the different b hadrons are varied in the Monte 



Carlo by their errors according to . 



c fragmentation: cc Monte Carlo events are reweighted by varying the Peterson frag- 
mentation parameter ep in the range of 0.025 < ep < 0.031 ||12|| . This corresponds to a 
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variation of the mean scaled energy 2(E c )/^/s of weakly decaying c hadrons in Z° decays 
in the range of 2(E c )/- s /s = 0.484 ± 0.008. 

c decay multiplicity: The average charged track multiplicities of D + , D° and D+ decays 
are varied in the Monte Carlo within the ranges of the experimental measurements ||27|| . 

c hadron composition: The D° fraction is written as /(D°) = 1 — f(D + ) — f(D+) — 
f(cbaryon)- The last three parameters are varied independently by their errors according 
to [0 to evaluate the uncertainty on the charm efficiency 

c hadron lifetimes: Charmed hadron lifetimes are varied within their experimental 



errors according to [21 



5.1.4 Four-fermion Background 



Above the WW production threshold, the four-fermion background is largely dominated by 
W pairs. The uncertainty in the W pair production cross-section is taken into account and 
has been found to have a negligible systematic effect. The background from W pairs has the 
highest probability to be accepted in the tagged sample when one or both W bosons decay into 
a final state containing a charm quark. The systematic error on the W pair tagging efficiency 
is estimated by varying the charm physics modelling parameters as described above. The effect 
of the detector resolution is also taken into account, as described in Section 15.1.6. 



5.1.5 Monte Carlo Statistics 

Tagging efficiencies and charge identification probabilities are varied by the statistical error 
arising from the finite number of Monte Carlo simulated events. 



5.1.6 Track Reconstruction 

The effect of the detector resolution on the track parameters is estimated by degrading or 
improving the resolution of all tracks in the Monte Carlo simulation. This is done by applying 
a single multiplicative scale factor to the difference between the reconstructed and true track 
parameters. A ±10% variation is applied independently to the rip and rz track parameters. 

In addition, the matching efficiency for assigning measurement points in the silicon mi- 
crovertex detector to the tracks is varied by 1% in rip and 3% in rz. The systematic errors 
resulting from the individual variations are summed in quadrature. 



5.2 Systematic Errors on A^ B and Ap B with the Lepton and Slow 
Pion Tag 

Three different groups of systematic errors have been considered: those from detector effects, 
those related to physics models and external inputs used in the analysis, and effects introduced 
by the limited Monte Carlo statistics and the fitting procedure. A list of all systematic errors 



is given in Table ITT 
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Error source 


A^FB 


A^ B 


Track reconstruction 


0.088 


0.070 


Lepton and slow pion identification 


0.009 


0.013 


Input modelling 


0.045 


0.022 


Fragmentation modelling 


0.024 


0.033 


Semileptonic decay models 


0.007 


0.005 


Branching ratios 


0.008 


0.006 


B — B mixing 


0.003 


0.001 


Background asymmetries 


0.024 


0.066 


Final state QCD corrections 


0.008 


0.011 


Monte Carlo statistics 
Fitting procedure 


0.061 
0.009 


0.032 
0.003 


Total systematic error 


0.122 


0.111 



Table 11: A breakdown of the systematic errors is given for the asymmetry measurement 
with leptons and slow pions at 189 GeV. Similar systematic errors have been obtained for the 
measurements at the other centre-of-mass energies above the Z° peak. 



5.2.1 Detector Systematics 

The lepton and slow pion identification relies on the proper modelling of the detector response 
in the Monte Carlo simulation. 

• Track Reconstruction: The influence of the simulated detector resolution and matching 



efficiency for hits in the microvertex detector is estimated as described in Section |5.1.6 . 

Lepton Identification: The fractions of misidentified electrons and muons are taken 
from Monte Carlo simulation. The modelling of the input variables to the artificial neural 
network for electron identification has been studied on Z° calibration data in a manner 



similar to that described in p7| . Differences between the data and the modelling in the 
Monte Carlo simulation have been determined using a pure sample of electrons from 
photon conversions and an inclusive sample of tracks depleted in conversion electrons. 
The dependence of the efficiency and background contamination of the selected electron 
sample on these differences has been studied, and the resulting uncertainties for each 
input have been added in quadrature. In addition, it has been shown that there is no 
large dependence of the resulting systematic error on the position of the A4i cut. The 
systematic error has been evaluated separately for each year of data taking at energies 
above the Z° peak; the uncertainty on the efficiency is around 10%, and that on the 
background contamination around 20%. 

The performance of the muon tagging network has been studied on 77 — ► events 
recorded at LEP energies above the Z° peak and on an inclusive sample of tracks from 
the Z° calibration data that fail the "best match" preselection criterion. The muon and 
background rates have been compared between data and Monte Carlo simulation in bins 
of Mjj,, and the simulation has been reweighted to match the data distribution. The 
reweighting factors in each bin have been varied independently by their statistical errors, 
and the resulting variations of the measured asymmetries added in quadrature. For the 
inclusive muon sample selected with a cut at A/" M > 0.65, the reweighting factors for muon 
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signal and background are 0.95 ± 0.05 and 0.97 ± 0.04, respectively. 

Conversion candidates are explicitly removed from the sample of events used in the fit. 
The M CY output distribution is well described in the simulation. The conversion finding 
efficiency is tested in the data using an algorithm which does not need the M e \ outputs 
for the two tracks as input. The accuracy of this test is 18%, and the conversion rate is 
varied by this amount to calculate the systematic error. The rate of non-prompt muons 



from pion and kaon decays in flight has been studied previously in P0[ , and has been 
found to be modelled to within 9%. 

The four above systematic errors are added in quadrature to yield the value given in 



Table 11 



Slow Pion Efficiency: Because of the large backgrounds in the slow pion sample, it 
is crucial to check that the slow pion reconstruction is modelled correctly in the Monte 
Carlo simulation. The efficiency has been studied by reconstructing slow pions in the 
Z° calibration data using the same cuts as described in Section [4.2.3| , except for the 



upper momentum cut which has been lowered to p™ x (\fs = rrizo) = 4.17 GeV. The 
slow pion content in this sample has been measured using a fit to the pi spectrum. 
Functions which are found to describe the shape of the background distribution well in 
the simulation are fitted to the pi sideband from 0.03 to 0.10 GeV 2 . The number of 
signal slow pion events is then determined by extrapolating the background estimate to 
Pt=0 GeV 2 . In the Z° calibration data, this fit has a relative statistical precision of 6%, 
but an additional systematic error of 23% is assigned to cover biases resulting from the 
extrapolation procedure and the particular choice of fit function. 

Combining this result with the OPAL measured yield of D* + mesons in hadronic Z° 
decays of n z o^ D *+x = 0.1854 ± 0.0041 ± 0.0059 ± 0.0069 g|, the efficiency of the slow 
pion reconstruction in Z° decays is measured to be (37.3 ± 9.0)%, where the relative 
errors of 6% and 23% have been added in quadrature. This is consistent with the value 
in the simulation of (31.9 ± 1.5(stat.))%. The fractional error on the slow pion efficiency 
is then assigned to the slow pion efficiencies at energies above the Z° peak, which are 
taken from the simulation. In order to determine the systematic error on the asymmetry 
measurements, the expected number of slow pions in each of the different subsamples is 
varied according to the error on the efficiency, and balanced by the expected number of 
background events so that the total number of tagged events is kept constant. 

No additional systematic error is assigned for the modelling of the shape of the slow 
pion signal pi distribution since the signal shape in Z° decays has been found to be well 
described in the simulation. 

Modelling of Artificial Neural Network Inputs: Each of the input distributions 
used in the flavour separation networks has been compared between data and Monte 
Carlo simulation. The simulated distributions are reweighted for each input variable in 
turn to agree with the corresponding data distributions, and the analysis is repeated with 
the weighted events. The observed differences from the original fit result are added in 
quadrature to yield the systematic uncertainty due to the modelling of the input variables. 
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5.2.2 Physics Systematics 



The fragmentation of charm and bottom quarks and the momentum spectra of leptons emitted 
in their semileptonic decay are described using phenomenological models tuned to experimental 
data. Systematic errors introduced by models of heavy hadron semileptonic decays, by the use 
of externally measured inputs, and from the assumed asymmetries of background candidates 
are assessed as follows: 

• Lifetimes: The lifetimes of weakly decaying b- and c-flavoured hadrons are varied as 



described in Section |5.1.3|. The resulting changes in measured asymmetries are negligible. 



Fragmentation: Bottom and charm fragmentation uncertainties are estimated as de- 



scribed in Section 6.1.3 



Semileptonic Decay Models: Systematic effects due to the modelling of semileptonic 



decays of heavy hadrons are studied following the recommendations in [12|. The lepton 
momentum spectra in the Monte Carlo simulation are reweighted to different theoretical 
models, with ranges of parameters chosen such that the experimental errors are covered. 

Branching Ratios and Hadronisation Fractions: The values for the semileptonic 
branching ratios B(b — > £), B(b — > c — > £), B(b — >c— >£), and B(c — > £) are taken 



from |L2|], and they are varied within their errors. Similarly, the values for the hadroni- 
sation fractions f(b — > c — > D* + ) and /(c — > D* + ) are taken from [22| and varied within 
their errors. 

In the absence of any measurements, the hadronisation fraction /(b — > c — > D*~) is taken 
from the simulation and varied by 100% in order to assess the systematic error. 

B — B Mixing: The B — B mixing parameter is taken as an external input. The value 



of x — 0.118 ± 0.006 [gl[] is used, and the variation within its error yields the systematic 
uncertainty listed in Table |TT| . 

Background Asymmetries: Prompt leptons and slow pions from four-fermion back- 
ground in the selected event sample can affect the measured asymmetries. They may 
lead to a rather large observed forward-backward asymmetry, but their contribution to 
the overall sample is small. The asymmetry from these events is varied within ±0.5 to 
determine the systematic error. The error due to the uncertainties on the bottom and 
charm asymmetries in radiative events is negligible. 

The asymmetry in the y distribution of mis-identified lepton and slow pion candidates is 
assumed to be zero for any event type. For lepton candidates, this assumption is verified 
within an accuracy of 2% on an inclusive sample of tracks that pass the lepton momentum 
cuts and do not pass the cut on P sig . For slow pions, a sample of tracks that pass the 
cuts on the slow pion momentum and transverse momentum, but do not pass the cut on 
"P S ig, is used. The asymmetry of slow pion background is found to be zero with an error 
of 4%. The asymmetries from lepton and slow pion backgrounds are then varied within 
these limits to assess the associated systematic errors. 



All the above errors are added in quadrature to yield the uncertainty given in Table ITT 



Final State QCD Corrections: As described in Section |5.1.2| , half the QCD correction 
to Ap B and Ap B as computed using ZFITTER is assigned as systematic error. 
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5.2.3 Monte Carlo Statistics 



An error arises from the limited statistics in the Monte Carlo simulation that is used to predict 
the fractions fi for each subsample (see Section 3~3). This uncertainty has been evaluated by 
varying for each subsample in turn the contribution from each source by its statistical error, 
redoing the fit, and adding all observed differences in quadrature. 



5.2.4 Fitting Procedure 



The fitting procedure has been studied in a large number of simulated experiments with the 
same statistics as in the actual measurements. It has been found to be essentially unbiased for 
the large data samples recorded at 183 GeV and 189 GeV. For the datasets at 133 GeV, 161 GeV, 
and 172 GeV, however, the statistical error from the fit has been found to be underestimated 
by a factor of up to 1.16 (bottom) or 1.05 (charm asymmetry), and biases of up to 0.12 a 
and 0.06 a have been found for the fitted bottom and charm asymmetry, respectively, where in 
each case a denotes the statistical error. The statistical errors are adjusted, the measurements 
are corrected according to the observed biases, and the full bias is treated as an additional 
systematic error. 



5.2.5 Cross-Checks 



Consistent results are observed when the fit is repeated with the numbers of bins in At,, A/b c , 
A4, and \y\ varied independently between 2 and 4. In addition, all selection cuts have been 
varied, with no significant deviations in the results. Consistent results have been found when 
the slow pion tagged events are not used in the fit. In the 189 GeV data sample, the statistical 
error of the charm asymmetry increases from 0.18 to 0.21 without the information from slow 
pion tagged events, while the error on the bottom asymmetry stays the same. 



6 Conclusions 



Using data collected at centre-of-mass energies between 130 GeV and 189 GeV with the OPAL 
detector at LEP, the relative e + e~ — > bb production rate and the forward-backward asymme- 
tries in e + e~ — > bb and e + e~ — > cc production have been measured to be 



Energy 




^FB 


^FB 


133 GeV 
161 GeV 
172 GeV 
183 GeV 
189 GeV 


0.190 ±0.023 ±0.007 
0.195 ±0.035 ±0.007 
0.091 ±0.034 ±0.005 
0.213 ±0.020 ±0.009 
0.158 ±0.012 ±0.007 


0.19±0.30±0.12 
-0.03 toil ±0.11 
0.82 ±0.14 
0.77 toll ±0.10 
0.63 t°oil ±0.10 


0.50 lg;l±0.13 
0.87 ±0.12 
0.69 ^±0.12 
0.55 1^±0.11 

o.5o ±8:ii ±0.11 



where in each case, the first error is statistical and the second systematic. These values are 
illustrated in Figure § together with the dependences on the centre-of-mass energy as predicted 
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in the Standard Model. For all measurements, good agreement is observed with the Standard 
Model expectation. The measurements of bb production presented in this paper supersede the 
previously published values of references [f§,§|]. 
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